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FINAL  TECHNICAL  REPORT  FOR  AFOSR  CONTRACT  F44620-76-C-0044, 

PROJECT  NO.  2307,  FOR  THE  PERIOD 
1  DECEMBER  1975  THROUGH  31  MAY  1980 

Comprehensive  List  of  Research  Objectives 

a.  Study  hypersonic  and  transonic  theory  with  emphasis  on  wing-body 
interference  and  other  nonplanar  problems. 

b.  Study  refraction,  reflection  and  diffraction  processes  occurring 
in  nonlinear  wave  propagation  in  steady  and  unsteady  flows. 

c.  Investigate  transonic  wind  tunnel  wall  interference  in  steady 
and  unsteady  flows. 

d.  Apply  inverse  solvers  to  slender  configurations  in  transonic 
flow. 

Substantive  Statement  of  Significant  Accomplishments  and  Progress  Toward 
Achieving  Research  Objectives 

Regarding  Items  a  and  b,  new  analytical  hypersonic  flow  solutions 
over  glide  vehicles  consisting  of  delta  wing-bodies  were  obtained  giving 
Insight  Into  the  important  role  of  such  geometric  features  as  the 
longitudinal  area  progression  of  the  body.  A  generalized  area  rule  was 
developed  in  which  the  change  in  L/D  of  the  clean  delta  wing  due  to 
body  addition  could  be  quantified  solely  in  terms  of  this  progression 
and  the  body  cross-sectional  shape.  Details  of  this  rule  are  given  in 
C3  in  the  Cumulative  Chronological  List  of  Publications  provided  in 
what  follows.  Further  information  with  emphasis  on  Item  b  will  be 
presented  in  C8.  These  publications  deal  with  the  nonlinear  wave  inter¬ 
action  processes  occurring  in  hypersonic  wing-body  interference.  Their 
exploitation  to  achieve  favorable  L/D  and  use  of  transparent  analytical 
relationships  derived  in  this  contract  is  described  in  C3. 

In  what  follows,  progress  related  to  Items  c  and  d  will  be  detailed. 


1 


Rockwell  International 

Science  Cantor 

SC5054.5FR 


Steady  Transonic  Tunnel  Wall  Interactions  and  Far  Fields 

A  Green's  function  method  has  been  developed  to  predict  the 
asymptotic  flow  upstream  and  downstream  of  airfoils  and  wings  In  wind 
tunnels.  Originally,  the  method  was  developed  in  a  Karman-Guderley  (KG) 
framework  and  in  more  recent  phases  has  been  extended  to  handle  modified 
small  disturbance  theory  (MSD)  flows.  The  analytical  results  of  these 
studies  have  been  reported  in  C4  and  C6 ,  and  further  details  and  numerical 
results  will  be  given  in  C9  and  CIO,  as  well  as  other  publications. 

As  compared  to  the  KG  formulation,  the  effect  of  the  MSD  modifications 
was  demonstrated  to  substantially  change  the  structure  of  the  volume  source 
term  appearing  as  the  dominant  term  in  the  far  field  for  the  perturbation 
velocity  potential.  Additional  contributions  were  demonstrated  to  occur 
due  to  the  influence  of  vorticlty  on  the  wing  and  wake  and  sources  over 
the  wing.  Analysis  of  the  appropriate  mathematical  representations  for 
these  effects  gave  rise  to  the  following  findings: 

1.  For  closed  trailing  edge  wings  the  pressure  associated  with 
the  wing  confined  between  solid  and  free-jet  walls  damps 
exponentially  far  upstream  and  downstream  of  the  wing. 

2.  For  open  trailing  edges  such  as  those  associated  with  jet 

flap  implementations  or  intermediate  inverse  (design)  solutions, 
a  source  term  evolves  with  solid  walls  which  is  proportional 
to  tue  Integral  of  the  trailing  edge  thickness  over  the 
wing  span  and  inversely  proportional  to  the  Glauert  compressi¬ 
bility  factor. 

3.  As  in  incompressible  flow,  for  large  distances  downstream, 
the  effect  of  the  shed  trailing  vortices  has  been  shown  to 

be  identical  to  that  produced  by  a  collection  of  two-dimensional 
vortices  of  individual  intensity  proportional  to  the  rate  of 
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change  of  the  spanwlse  load  distribution.  The  location  of 
each  of  the  two-dimensional  vortices  is  on  the  projection  of 
the  corresponding  trailing  vortex  in  the  Trefftz  plane.  For 
rectangular  tunnel  cross-sections,  the  appropriate  expression 
for  the  induction  is  a  double  Fourier  series  in  the  transverse 
coordinates  normalized  by  the  characteristic  height  and  breadth 
of  the  tunnel  dimensions.  This  can  be  also  expressed  in  terms 
of  elliptic  theta  functions.  The  effect  of  the  transonic 
nonlinearity  is  exhibited  in  the  higher  order  terms  and  the 
near  field  coupling  with  spanwlse  loading.  This  validation 
of  the  transonic  application  of  Blot  Savart's  law  can  be  used 
to  readily  estimate  the  asymptotic  upstream  and  downstream 
behavior  associated  with  tunnel  cross  sectional  shapes  other 
than  rectangular  for  free-jet,  slotted,  and  solid  wall  boundaries. 
This  can  be  achieved  through  the  use  of  mapping  and  Image 
reflection  principles.  For  porous  walls  the  procedure  reverts 
back  to  the  Green's  theorem  method.  Here,  an  unconventional 
Sturm  Liouvllle  problem  arises  in  which  the  eigenvalues  appear 
in  the  boundary  conditions. 

4.  For  closed  trailing  edges,  the  exponential  damping  exhibited 
by  the  tunnel-wall  far  field  solution  contrasts  with  the 
algebraic  behavior  associated  with  the  free-fleld  configuration. 
For  open  trailing  edges,  the  associated  source  effect  damps 
algebraically  with  distance  in  contrast  to  the  growth  behavior 
exhibited  for  the  confined  wall  case  Indicated  in  (2). 

More  recent  phases  of  the  effort  deal  with  assessment  of  tunnel  wall 
Interference  on  compound  lifting  units.  Typical  results  from  C7  are 
shown  in  Fig.  1.  These  show  the  effects  of  solid  walls  on  the  mutual 
aerodynamic  interaction  between  dual  lifting  elements  representative  of 


3 


<!» 


Rockwell  International 

Science  Center 

SC5054.5FR 


—  FREE  FLIGHT 
—SOLID  WALL  nQ 


OC-0.4  + 

<  Ji 

3-02* 
o  ! 

i 


|  Cy 


Cw— 4 
— s — J 


Rockwell  International 

Science  Center 


SC5054 . 5FR 


a  two-dimensional  characterization  of  canard-wing  interference.  These 
calculations  used  special  far  fields  developed  for  multiple  systems. 

The  chordwise  pressure  distributions  indicate  that  the  wall  choking 
effect  increases  the  rear  lifting  surface  (wing)  effective  incidence 
angle,  apparently  resulting  in  a  larger  supercritical  region  and  the 
addition  of  a  terminating  shock.  On  the  other  hand,  the  forward  lifting 
surface  (canard)  experiences  a  lower  incidence  angle  compared  with  the 
free  flight  case,  and  the  shock  is  much  weaker.  These  striking  trends 
could  reverse  with  other  combinations  of  the  parameters  such  as  gap, 
stagger,  Mach  number  and  incidence  but  are  felt  indicative  of  variations 
of  an  important  subclass  of  these  variables.  Additional  far  field 
analytical  results  which  apply  to  slotted  and  free  jet  cases  are 
described  in  C7  with  a  summary  of  the  analysis.  Details  will  be 
given  in  C9  and  CIO.  Some  of  these  results  have  been  presented  at 
the  recent  Adaptive  Wall  Symposium  at  A£DC  in  17  and  relate  to  important 
difficulties. in  accounting  for  multi-wing  interactions  in  the  convergence 
of  the  adaptive  wall  iteration  process  in  obtaining  interference  free 
results  for  the  tail  chordwise  pressures. 

A  special  phase  of  the  research  effort  has  been  devoted  to  the 
nonself-ad joint  problems  associated  with  modeling  porous  walls.  These 
studies  which  will  be  reported  in  CIO  are  being  used  to  assess  blockage 
effects  and  can  be  generalized  to  the  treatment  of  adaptive  tunnels. 
Because  of  a  difficulty  associated  with  an  eigenvalue  appearing  in  the 
boundary  conditions  for  the  porous  case,  special  procedures  were  developed 
in  FY  1979  to  determine  the  Green's  function  which  appears  as  a  kernel 
in  the  integral  representation  for  the  perturbation  potential  and  is  a 
solution  to  the  nonself-adjoint  boundary  value  problem.  Our  new  method 
utilizes  auxiliary  functions  which  form  a  "dual  basis"  of  eigenfunctions 
having  self-adjoint  properties  and  orthogonality.  With  this  procedure, 
the  Green’s  function  has  been  developed  for  circular  test  sections  as 
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an  example  of  the  approach  which  has  been  conceptualized  for  arbitrary 
cross  sections.  This  work  will  serve  as  a  basis  for  the  continuation 
of  the  effort  dealing  with  a  systematic  asymptotic  procedure  for  the 
assessment  of  tunnel  wall  blockage  effects,  expanding  the  focus  to  near 
field  aspects  in  addition  to  the  downstream  and  upstream  behavior 
currently  being  investigated. 

Another  portion  of  the  effort  has  involved  sensitivity  analyses 
regarding  the  near  field  impact  of  errors  in  the  satisfaction  of  the 
correct  far  field  boundary  conditions  in  a  wind  tunnel  environment. 

This  activity  is  partially  motivated  by  adaptive  wall  applications. 

For  the  linear  case,  a  perturbation  formulation  has  been  used  to  show 
that  an  error  in  the  estimate  of  the  far  field  produces  the  same  order 
of  effect  in  the  near  field.  For  the  nonlinear  case,  this  generalization 
is  not  immediately  evident,  and  numerical  methods  are  currently  being 
utilized  to  provide  information  regarding  the  effects  of  such  parameters 
as  wall  height,  chord  and  control  section  lengths. 

Typical  results  from  these  two-dimensional  studies  are  shown  in 
Figs.  2-4.  For  these  calculations,  a  NACA  0012  was  modeled  in  a  solid 
wall  tunnel  using  a  two-dimensional  SLOR  solver  for  the  Karman-Guderley 
transonic  small  disturbance  model.  The  far  field  was  computed  using 
expressions  from  C7.  In  Fig.  4,  confined  lift  coefficients  are  plotted 
as  a  function  of  the  streamwise  location  of  the  upstream  and  downstream 
control  boundaries  situated  at  x„  chord  lengths  from  the  center  of  the 

f 

airfoil  which  is  situated  on  the  tunnel  centerline.  The  calculations, 
which  were  performed  at  ■  0.75  and  a  ■  2®,  retaining  only  the  dominant 
terms  in  the  far  field  indicate  that  the  appropriate  Xp  «  00  asymptotic 
value  is  approached  more  rapidly  for  lower  height  to  chord  ratios,  H. 

This  trend  is  symptomatic  of  the  transition  of  decay  of  disturbances 
between  small  and  large  H.  For  smaller  values  of  this  parameter,  our 
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far  field  research  has  indicated  that  the  decay  is  exponential,  whereas 
for  large  values  it  becomes  algebraic,  appropriate  to  the  free  field. 

This  behavior  will  be  a  crucial  part  of  the  matching  aspects  in  the 
asymptotic  tunnel  interference  theory. 

In  Fig.  3,  similar  studies  were  performed  indicating  the  approach 
to  the  asymptotic  values  of  the  confined  lift  coefficients  with  a 
second  order  expression  for  the  far  field.  Here,  it  is  evident  that 
inclusion  of  more  terms  speeds  up  the  approach  to  the  asymptotic  C^. 
Accordingly,  the  dimensions  of  the  computational  domain  can  be  reduced. 

For  adaptive  walls,  the  significance  of  this  result  is  to  reduce  the 
streamwise  extent  of  the  control  domain.  Figure  4  reaffirms  these 
trends  at  a  »  4*  and  also  Indicates  the  larger  x  required  to  reach 
the  asymptotic  values.  This  trend  is  associated  with  the  stronger 
disturbances  produced  by  the  larger  incidence. 

To  study  analogous  trends  in  three-dimensional  flows,  parametric 
studies  were  performed  using  the  geometry  Indicated  in  Fig.  5.  Here, 
a  rectangular  planform  wing  with  a  NACA  65A006  section  was  simulated  in 
a  rectangular  cross  section  solid  wall  tunnel.  Calculations  were  performed 
for  varying  M^,  a,  H,  D,  and  wing  aspect  ratio,  studying  the  decay 
of  the  pressure  field  with  x^  in  this  investigation.  For  the  geometry 
considered,  the  far  field  analytic  formula  further  developed  in  recent 
phases  of  the  current  contract  can  be  written  in  the  form  (using  the 
notation  of  Fig.  5) 

♦  *  ^FF  “  *W  +  Xwa  +  TT  +  Xvs  as  *  ±m  (la) 


where 
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■  (a)'  *  (*)' 

"  (1  +  5n0)(1  +  60m)  *  n  +  m  *  0 

-  Kronecker  delta  function  . 

Equations  (11)  can  be  sunmarlzed  In  the  form: 

<J>FF  =  A(y,z)H(x)  +  Be  sin  Xy  as  |x|  ■+  “  (2) 

where  A  and  B  can  be  obtained  by  comparison  with  (11) . 

Some  typical  results  of  the  parametric  studies  will  now  be  discussed. 
In  Fig.  6,  calculations  for  the  downstream  decay  of  the  absolute  value  of 
pressure  on  the  bottom  wall  |$x(x,-H,0) 1  are  shown  with  a  simplified 
version  of  (12),  <t>x(x,y,±20c)  ■  0,  imposed  at  the  computational  upstream 
and  downstream  boundaries.  For  the  cases  to  be  discussed  hereinafter, 

■  0.9  and  a  »  5°.  The  correct  decay  rate  as  shown  by  the  dash  dot 
line  from  (12)  is  compared  to  the  numerically  obtained  result,  indicating 
that  placing  the  infinity  boundary  condition  at  a  finite  upstream  and 
downstream  location  still  produces  exponential  decay  but  with  a  decay 
rate  different  than  X. 

Figure  7  indicates  the  effect  of  utilizing  the  two  term  asymptotic 
Eqs.  (1)  in  the  far  field.  For  an  intermediate  range  of  x,  it  is  evident 
from  the  plot  that  the  wall  pressure  distribution  correctly  tracks  the 
proper  decay  rate  for  xp  sufficiently  large,  l.e.,  20.  The  lack  of 
correlation  near  the  wing  is  expected  since  that  region  is  nonlinear. 

This  behavior  gives  an  indication  of  the  control  length  limits  for  adaptive 
wall  applications.  It  is  furthermore  significant  that  anamolous  results 
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Upstream  and  downstream  decay  of  absolute  value 
of  pressure,  H  *  D  ■  4c,  Mm  ■  0.9,  a  ■  5*, 

NACA  65A006  section,  AR  -  2,  $  (x,y,±20c)  -  0. 
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are  obtained  for  x  large  where  a  new  phenomenon  takes  place.  Here, 
discretization  errors  resulting  from  the  expanded  grid  overshadow  the 
alleged  improved  accuracy  of  the  Prandtl-Glauert  asymptotic  (1) .  A 
more  correct  picture  than  that  obtained  from  the  analytical  far  field 
analysis  leading  to  (1)  would  evolve  from  a  study  of  the  relevant 
difference  equation  properties  rather  than  the  continuum  limit  as 
x  ■*  ®.  This  limit  is  extremely  subtle  involving  a  combination  of  the 
expanding  mesh  and  large  x.  As  for  the  continuum  formulation,  the 
discrete  problem  is  linear.  Appropriate  solution  procedures  are 
currently  under  study. 

In  Figs.  8  and  9,  the  significance  of  the  far  field  in  regard 
to  the  value  of  the  lift  and  drag  coefficient  on  the  rectangular  wing 
for  the  previous  values  of  the  parameters  is  shown.  Here,  the  free 
field  coefficient  is  denoted  on  the  ordinate  scale  by  the  indicated 
arrows.  From  these  results,  it  is  obvious  that  the  position  of  the 
control  boundary  has  a  decisive  effect  on  the  magnitude  of  the  force 
coefficients.  Studies  such  as  these  give  rise  to  the  striking  conclusion 
that  the  location  of  the  far  field  downstream  control  boundary  can  affect 
the  lift  coefficient  by  an  amount  exceeding  the  basic  solid  wall  blockage 
correction  to  the  free  field  values.  This  emphasizes  the  significance  of 
the  proper  application  of  downstream  and  upstream  far  field  conditions 
in  adaptive  and  nonadaptive  tunnel  Interference  modeling. 

Unsteady  Tunnel  Wall  Interference  Effects 

For  unsteady  wind  tunnel  flows,  a  mathematical  framework  schematically 
Indicated  in  Figs.  10  and  11  has  been  devised  to  provide  a  consistent 
approximation  procedure  for  the  characterization  of  the  far  fields.  This 
approach  corresponds  to  the  unsteady  counterpart  of  the  Integral  repre¬ 
sentation  method  utilized  for  steady  confined  flows.  Results  obtained 
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Fig.  LO.  Far  field  for  dilating  profile. 
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.V.3.:  FOR  w-0,  NEED  ONLY  A  FEW  TERMS  IN  SUMS  -  OTHERWISE  FAR  FIELDS 
ARE  PER 1001 C  ANO  LARGE. 

Fig.  11.  Dilating  airfoil  in  presence  of  solid  walls. 
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in  FY  1979  are  for  harmonically  dilating  transonic  airfoils  between 
solid  walls  as  well  as  in  a  free  field.  The  structure  obtained  confirms 
the  lack  of  damping  of  the  disturbances  for  moderate  frequencies  as 
compared  to  the  steady  case.  A  similar  result  was  obtained  in  FY  1978 
using  a  heuristic  model  of  a  lifting  airfoil  as  a  pulsating  vortex. 

Extensions  of  the  systematic  approximation  procedure  to  treat  lifting 
airfoils  and  wings  using  Kirchoff's  formula  and  the  present  method  are 
in  progress. 

Studies  of  Transonic  Slender  Body  Flows 

An  inverse  method  is  being  applied  to  design  slender  bodies  to 
support  preassigned  pressure  distributions  with  the  eventual  application 
to  cruise  missiles  and  slender  airplane  shapes  such  as  fighter  configurations. 
In  Fig.  12,  results  are  shown  for  an  axisymmetric  solver  which  is  being 
applied  in  the  design  mode.  Special  boundary  discretization  methods 
have  been  devised  to  provide  consistency  in  the  inverse  mode.  The 
excellent  agreement  with  the  experimental  results  of  C8  attest  to  the 
accuracy  of  our  procedures.  Currently  an  iterative  algorithm  is  being 
checked  out  for  the  inverse  problem.  An  outline  of  the  procedure  is 
given  in  Fig.  13.  Here,  R(x)  is  the  body  radius,  g'  is  the  nonlinear 
function  related  to  the  surface  pressure  coefficient,  T  is  the  fineness 
ratio,  and  n  is  an  iteration  counter.  The  slender  body  theory  employed 
for  this  phase  has  been  generalized  to  provide  indications  of  the 
existence  of  an  unsteady  equivalence  rule  which  represents  a  new 
generalization  of  the  steady  theorem  devised  by  Oswatitsch.  Besides 
pointing  to  unconventional  beneficial  interfering  concepts  for  reduction 
in  wave  drag,  the  results  of  this  analysis  could  find  application  in 
cost  effective  means  of  calculation  of  dynamic  longitudinal  stability 
properties  such  as  damping  in  pitch  of  slender  vehicles  and  missiles  in 
the  transonic  regime.  An  important  cost  savings  will  be  achieved  by  the 
equivalence  rule's  removal  of  one  of  the  spatial  variables  of  the  problem. 
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12.  Comparison  with  experiment  of  present  slender  body 

analysis  code  with  design  oriented  boundary  discretization.  . 
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Fig.  13.  Flow  chart  for  slender  body 
iteration  procedure. 
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Aerodynamics  Simulation  Facility.  This  committee  has  Dr.  W.F.  Ballhaus  \ 

as  its  chairman  and  meets  approximately  twice  yearly.  Other  individuals  [ 

on  this  committee  include  P.  Rubbert,  R.  Melnik,  W.  Hankey,  and  F.R.  Bailey. 

■i 

Other  Interactions 

Dr.  Malmuth  is  chairman  of  the  Rockwell  fluid  dynamics  technical 
panel  and  serves  on  the  AIAA  Fluid  Dynamics  Committee. 

Other  Statements  | 

The  thrust  of  the  research  effort  is  to  couple  asymptotic  methods  I 

with  computational  techniques  to  enhance  the  efficiency  of  the  latter  | 

in  analyzing  nonlinear  flows.  Results  obtained  from  the  far  field  g 

activity  and  application  of  new  nonlinear  slender  body  design  methods  I 

to  be  developed  in  the  next  phase  of  this  contract  should  provide  the  | 

basis  of  handling  three-dimensional  configurations  to  achieve  marked  | 

reductions  in  the  cost  of  transonic  aerodynamic  design  and  analysis.  I 

In  this  regard,  the  consistent  inverse  design  procedure  developed  under  | 

Rockwell  IR&D  and  reported  in  C4  and  C6  will  be  applied  in  this  effort  I 

to  achieve  optimum  configurations.  An  immediate  fallout  of  this  effort 
should  be  applications  to  cruise  missile  design.  Other  derivatives 
of  the  work  will  be  the  development  of  technology  baselines  for  the 
design  of  the  next  generation  of  energy  efficient  and  high  performance 
tactical  and  strategic  vehicles. 
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